Band Structure and Accumulation Point in the Spectrum of Quantum Collision 
Operator in a One-Dimensional Molecular Chain 
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We consider the eigenvalue problem of a kinetic collision operator for a quantum Brownian par- 
ticle interacting with a one-dimensional chain. The quantum nature of the system gives rise to a 
difference operator. For the one-dimensional case, the momentum space separates into infinite sets of 
disjoint subspaces dynamically independent of one another. The eigenvalue problem of the collision 
operator is solved with the continued fraction method. The spectrum is non-negative, possesses 
an accumulation point and exhibits a band structure. We also construct the eigenvectors of the 
collision operator and establish their completeness and orthogonality relations in each momentum 
subspaces. 
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I. INTRODUCTION 

Since the seminal work of Boltzmann's on dilute gas 
[1], the kinetic equation approach had been extensively 
applied to various fields, for example, the transport phe- 
nomena in solid state physics [2, 3], dense fluids and 
plasma [4, 5], optical transition in quantum optics and 
quantum information [6, 7], the motion of quantum 
Brownian particle in a potential [8], and etc. The ki- 
netic equation incorporates the effect of fluctuation on 
a test particle by considering the collision between the 
particle with its environment. 

A density matrix formulation for the relaxation process 
is desired since dissipation inevitably brings a pure state 
into a mixed state, which is incompatible with the unitary 
time evolution on the wave function level through the 
Schrodinger equation. This results in a non-hermitian 
collision operator that governs the relaxation of the den- 
sity matrix of the particle towards equilibrium. It is 
therefore important to study the eigenvalue problem of 
the collision operator to understand the relaxation pro- 
cess. 

As has been shown in Ref.[9], the eigenvalue prob- 
lem of the collision operator is closely related with the 
eigenvalue problem of the Liouville operator. However, 
spectral properties of the Liouville operator, as well as 
spectral properties of the collision operator, are not yet 
well-analyzed as compared with the spectral properties 
of the Hamiltonian. Indeed, despite the versatility of the 
kinetic equation approach and the existence of various 
methods in studying it [4, 5, 8], only a few models are 
known to be exactly solvable, such as, e.g., the Maxwell 
molecules [4] and the particle- field linear coupling model 
in one dimension (ID) [10]. The model that we are going 
to study is another example of an exactly soluble model 
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in ID on the level of kinetic equation. It gives rise to a 
few interesting results to be discussed below. 

We study a quantum Brownian particle interacting 
with a ID chain. The model is derivable, for instance, 
from an electron interacting with a background ID ion 
lattice through a deformation potential [11, 12], or a 
Davydov soliton produced by peptide oscillation prop- 
agating under the influence of a ID molecular chain [13]. 
Previous studies on this model were mainly carried out 
on the wave function level [13]. Here, we focus on the re- 
laxation process on the level of momentum distribution 
function of the particle under the thermal fluctuation of 
the background ID chain. 

The collision operator of this system exhibits a few in- 
teresting features [14]. In our quantum system, it takes 
a form of a difference operator, in sharp contrast to the 
usual differential operator in classical systems, as well as 
linear coupling quantum systems. Moreover, the dissipa- 
tive effect introduced by the collision operator is a purely 
quantum effect similarly to the case of the kinetic equa- 
tion in ID classical gas. Our collision operator possesses 
also a special character owing to its ID nature, i.e., the 
momentum states are coupled discretely, resulting in in- 
finite sets of dynamically disjoint momentum subspaces. 

The difference equation in the momentum space ren- 
ders the usual method in solving differential equations 
(such as eigenfunction expansion and the Green's func- 
tion method) inapplicable to this case. We instead find 
that the eigenvalue problem is best solved with the con- 
tinued fraction method [8, 15], which can be efficiently 
implemented numerically. 

The main results of the paper are as follows. We find 
that the spectrum of the difference operator is discrete in 
each disjoint momentum subspace and non-negative. It 
possesses a maximum eigenvalue in each momentum sub- 
spaces and an accumulation point. Moreover, the eigen- 
values are non-degenerate except in the steady mode, 
thus giving rise to a band structure. 

In Ref. [16], we presented some of the main features of 
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the solution by continued fraction method and the spec- 
trum, but without much discussion on the mathematical 
details. Here, we clarify the mathematical structure of 
the eigenvalue problem of the collision operator. While 
Ref. [16] used a symmetrized version of the collision op- 
erator, here we work with the original non-hermitian col- 
lision operator. We clarify the biorthogonal space struc- 
ture [15] of the operator that belongs to the generalized 
space [17]. We also introduce a scalar product to the 
Liouville space, clarify the detailed structure of the spec- 
trum and construct the eigenvectors. 

The discussion is organized as follows. In Section II 
we introduce the model and the difference collision op- 
erator. The disconnectedness nature of the momentum 
space is also explained. In Section III wc consider the 
eigenvalue problem and clarify the biorthogonal structure 
of the eigenvector. We then use the continued fraction 
method to solve the eigenvalue problem in Section IV and 
discuss the details of the spectrum in Section V. In Sec- 
tion VI we construct the cigcnhmctions of the collision 
operator and make a comparison with the approximate 
solution obtained by truncating the difference collision 
operator to be a second order differential operator. In 
Section VII we present conclusions. 



a periodic boundary condition leading to discrete mo- 
mentum p and wave numbers q with p/h,q = 27rj/L, 
where j = 0,±1,±2,---, respectively. 

The time evolution of the total system follows the 
Liouville-von Neumann equation. 



d 

with the Liouvillian £ defined by 



(3) 



(4) 



where p{t) is the density matrix of the total system, and 
the Liouvillian C is defined \>y Cp = [H,p\/h. We are 
interested in the time evolution of the reduced density 
matrix of the particle, and we average over the phonon 
field, 



f{t) = Trphp(t). 



(5) 



We assume that the phonon system is in thermal equi- 
librium with a temperature T represented by 



II. DIFFERENCE COLLISION OPERATOR 



(6) 



A. Model 

We consider the relaxation dynamics of a quantum 
Brownian particle interacting with the acoustic phonon 
field of a ID lattice [11, 12], 

P Q 

+ -/t E ^> + IK + o-l,) , (1) 

where wc restrict our consideration to the one particle 
sector of the Brownian particle. The momentum state 
vector of the particle is labeled by \p). The creation and 
annihilation operators of the phonon field with wave vec- 
tor q are a J and Ug, respectively. The length L is the 
dimension of the lattice. The particle has mass m and 
energy Sp = p^ /2m, whereas the acoustic phonons as- 
sume the dispersion relation ooq = c\q\. The particle is 
coupled to the phonon field through a deformation po- 
tential with efi^ective coupling [3, 18] 



Vn 



AwpMCOq ' 



(2) 



in which Aq is the deformation potential and pM is the 

mass density of the chain. Wc consider a weakly cou- 
pling case between the particle and phonons. We impose 



where /3 = l/ksT with the Boltzmann constant ks- 

We are interested in the limit L — >■ oo, where the wave 
number becomes continuous and the summation over the 
momentum is replaced by an integration. 



27r 



E 



dp, 



27r 



6{p-p'). (7) 



We will continue to use the discrete notation in our pre- 
sentation for compactness of expression. 

For the weak coupling system, the efi'ect of the inter- 
action between the particle and phonons can be approx- 
imated to the second order with respect to the potential 
Vq, and the time evolution equation of the momentum 
distribution function defined as f{p, t) = {p\f{t)\p) obeys 
a Markovian kinetic equation in a form idf{p,t)/dt = 
'iff{p,t), where "if is the generator of the time evolution 
(see e.g. Ref.[19]). 

In view of the fact that the eigenvalues of ^ are pure 
imaginary (which will be shown in due course), we will 
consider instead the operator K. = — z^, where K. is the 

collision operator that we will focus on from now on. 

The reduced dynamics follows a Markovian kinetic 
equation [14] 



|/(p,t)=^/(p,t), 



(8) 



where the action of the collision operator on the distri- 
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bution function is given by 



hq ^ p 



h 



X {[n{q) + l]f{p,t)-n{q)f{j>-hq,t) 



X [n{q)f{p, t) - [n{q) + + hq, t)^ | 



(9) 



The number density of the phonon obeys the Bose- 
Einstein distribution n(%) = l/[exp(/3fiwg) — 1]. We note 
that as h approaches zero, the coUision operator vanishes 
[14]. Therefore, the dissipation caused by the collision 
operator is a purely quantum effect. 

Notice that the collision operator takes the form of a 
difference operator. It can be expanded in an infinite 
series of higher derivative with respect to momentum, 
/(p ± hq) = exp{±hqd/dp)f{p). This is in contrast to 
the usual stochastic differential kinetic equation, such 
as the Fokkcr-Planck equation [8], which contains dif- 
ferential operators up to the second order in momentum. 
We cannot in general approximate the difference kinetic 
equation by truncating the series of differential operators 
to more than the second order, otherwise the positivity of 
the reduced distribution function cannot be maintained 
[8, 20]. 



B. Disjoint momentum subspaces 

From the resonance conditions represented by the delta 
functions in Eq. (9) for the absorption and emission of 
phonons, we obtain Sp' — Ep = hwg > and Sp' — Sp = 
— hojq < 0, respectively. Consequently, two adjacent mo- 
menta p and p' are related by 

- ± (p' _ p) = l{p'-p){p'+p±2) = 0, (10) 
where we introduce the dimensionless momentum 



p = p/{mc) , 



(11) 



and hq/{mc) = p' — p- Thus each p is coupled with two 
other momenta —p± 2. If we start with some po, then all 
the momenta coupled directly or indirectly to po can be 
reached by recursively applying the following formula. 



Pi±i = -Pi±2{-iy, 



0,±1,±2,.... 



(12) 



The solutions of this recursive formula are the mutually 
disjoint momentum subspaces each represented by po, 



Po;±i = i-iy{poT2i), i = 0,1,2,.... 



(13) 



Fig. 1 illustrates the disjoint momentum subspaces for a 

few values oi pQ. When pQ = 0, the positive and negative 
branches are mirror image of each other along the po = 




FIG. 1: Disjoint momentum subspaces, for (i) po = 0, (ii) 
— 1 < Po < 1 and (iii) po = 1. Solid lines connect the positive 
branch, and dotted lines connect the negative branch. 



axis, see Fig. l(i). At po = ±1) both branches will merge, 
see Fig. l(iii) for the case of po = 1- By varying po across 
the interval —1 < po < 1, all the momentum states are 
covered. The resonance conditions therefore divide the 
momentum space into infinite disjoint sets {po-i}, each 
set is represented by po in the interval — 1 < po ^ 1- 

The existence of the disjoint momentum subspaces is a 
special feature of a ID system. In higher dimensions, the 
energy and momentum conservation equations will give 
rise to, in the two dimensional case for instance, a set of 
closed lines. These lines may intersect with each other. 
As a result, all the momentum states become connected. 



III. EIGENVALUE PROBLEM OF COLLISION 
OPERATOR 

A. Biorthogonal basis 

By specializing to a particular momentum subspace 
{Po;i} connected to po, let us consider the eigenvalue 
problem 



K^4>^iiPo■,±i) = -A^(po)0^(po;±i) : 



(14) 



where we anticipate the fact that the eigenvalue is a func- 
tion of Po- We will label the components of the eigenvec- 
tors as 

(pt.iPo) ^ MPo;±i) , i = 0,l,2,.... (15) 
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When po lies in the interval (—1, 1), the eigenvalue prob- 
lem takes the form 

-A^</'p;0 = [I + ni2 + 2poM-i 

- [n(2 - 2po) - n{2 + 2po)]0M;O 

+ [1 + n{2 - 2poM+^ , \po\<l, (16a) 
= ^(2|Po;±.|-2)(/<±,_i 

- [1 + n(2|po;±^| - 2) + n(2|po;±i| + 2)]0±., 
+ [l + n(2|po;±.|+2)],^±.,+i, 

l<|PO;±i|, i = 1,2,3,..., 

(16b) 



where we define the dimensionless eigenvalue 



(17) 



and omit the po dependence on 0^. ■ to simplify the no- 
tation. The corresponding dimensionless time is then 
t = mAl/{pMh^c)t. 

For the boundary value po = ±1, Eq. (16a) will be 
modified. In this case, there is overlap in the momentimi 
states, pQ-ij-i =Po.i, see Fig. l(iii). As a result, (j^'^.i^-i = 
(/>~.j. Therefore, the components of the eigenvectors that 
need to bo considered reduce by almost half compared 
to the case po ^ ±1. At this point, we make use of the 
freedom to choose, as independent components of the 
eigenvectors, elements of the negative branch when 
Po = 1, and elements of the positive branch ^+ when 
Po = —1. We then have 

-A^^Mio = -?i(4)0m;o + [1 + ni'iMJ.i , Po = ±1 , 

(16a') 

whereas the rest of are still related to (p^-o via 
Eq. (16b). 

The eigenvalue problem can now be compactly written 
as a matrix equation, 

K{po) ■ 0^(Po) = -MPo)4>M) ■ (18) 
The eigenvector 0^ is a column matrix, 

fPnipo) = 4>Mpo) , (19) 



V 



whereas the collision operator 1<" is a tridiagonal ma- 
trix, with components Kij that can be read directly from 
Eqs. (16a) or (16a'), and (16b). They are 

Ko,o = -n{2 - 2po) - n(2 + 2po) , |po| < 1 , (20a) 
Ko^±x = 1 + n(2 T Po) , Ipo I < 1 , (20b) 

and 



K±i,±r = -1 - n(2|po;±.| - 2) - n(2|po;±»| + 2) , (20c) 

(20d) 
(20e) 



^±i,±(i+l) = l + «(2|pO;±z|+2). 

^±i,±(i-i) = n{2\pQ.±i\ - 2) , 



For Po = ±1, in place of Eqs. (20a) and (20b), we have 



ifo,o = -n(4) , 
ifo,±i = 1 + n(4) , 



Po = ±1 , 
Po = ±1 . 



(20a') 
(20b') 



The canonical equilibrium distribution function 
exp(— /9pQ._|_j/2) is the steady state of the collision op- 
erator with zero eigenvalue Ao = 0, where ^ = 1/T and 
T is the dimensionless temperature 



f = kBT/imc^) , 



(21) 



This can be inferred from Eqs. (A2) and (A3) in the 
Appendix A, with the help of relation Eq. (28) below. 
We denote the components of the steady state vector by 



-«,±./2 



With the normalization factor 



0,1,2,... . 



Z(po) = J2 



/2 



(22) 



(23) 



we have ip = <pQ. 

Due to the non-hermiticity nature of K, we also need 
to consider the left eigenvalue problem xj, • if = — A^xji 
[15], or equivalently, 

K\po) ■ xM) = -Mpo)Xi,{Po) ' (24) 

where is the left eigenvector of K. It can be shown 
that the components of K satisfy the relation (no sum- 
mation over i and j) 



In matrix form, we write this as 



(25) 



(26) 



where the o-product is defined by the right hand side 
(RHS) of Eq. (25). It follows from Eq. (18) that 

■ {cp-^ o <^^(po)) = A^(po)<^-^ o <^^(po) . (27) 
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A comparison with Eq. (24) shows that the left eigenvec- 
tors are related to their right counterparts by 

= o , (28) 

or equivalently, 

<P^. = 'P°X^■ (29) 

The left and right eigenvectors are biorthogonal [15]. 
This means that they are orthogonal and complete in the 
sense of 

xliPo) ■ MPo) = S^. , (30) 
E<^m(^o)-4(Po)=I, (31) 

respectively, where I is an infinite dimensional identity 
matrix. It should be noted that Eqs. (30) and (31) apply 
to each po subspace, and there is a separate complete set 
of orthonormal eigenvectors for each po subspace. 

In this space, we define the inner product of two vectors 
g and / as 

{{9{Po)\f{Po)))^['p-'og{Po)]^.f{Po) 

OO 

= E ^r'A*(Po)/.(Po). (32) 

i=—oo 

This inner product should be compared to the usual inner 
product in continuous space, 

/OO 
dp<f-^{p)g*{p)f{p). (33) 
-OO 

Eq. (30) is then equivalent to the norm of the right eigen- 
vectors, = S^y. Using the orthogonal relation 
(30), we show in Appendix A that the eigenvalues are 

non-negative. 

It is a common feature in dissipative systems that the 
eigenvectors belong to the generalized space[17], i.e., they 
are not objects in the Hilbert space. In other words, the 
pair of vectors arc generalized vectors. The 

vectors belong to the space of test function and is 
normalizable. On the other hand, belong to the dual 
space of the test function and arc not normalizable. For 
example, consider the steady mode that has components 
XO;i = 1 for all I. Xa is therefore non-normalizable under 
the inner product (32). 

To some extent, the non-hermiticity of the operator K 
is artificial, since one may transform it to a hermitian 
operator through a similarity transformation. A sym- 
metrized version of K as introduced in Ref. [14] achieves 
this aim. The transformed left and right eigenvectors in 
this case become normalizable. 



B. Time evolution 

We can determine the time evolution of an initially 
continuous distribution function f(p) by first discretizing 



it into a set of column vectors /(po) (see Eq. (19) for an 
explicit form of the vector), for all pQ. The components 
of the vector f{po) are labeled by 

/±(P0)=/(P0;±^), i = 0,l,2,--.. (34) 

Using the orthonormality of the basis vector (30), we 
can expand /(po) in terms of <p^{po). The time evolution 
of f{po) can then be worked out, 

/(po;t) = e^*-/(po;0) 

= EcM(Po)e-^''(^"'')V^(po), (35) 

where the expansion coefficient is 

c^{Po) = xliPo)- f{Po;0). (36) 
The time evolution of the inner product is 

{{f{Po;t}\f{Po;m = E |c^(Po)l'e-2^''(^"'')* , (37) 
ft 

after using Eqs. (28) and (30). It is a monotonously 
decaying function, in accordance to the Markovian nature 
of the kinetic equation under the A^t-approximation. 

Since the collision operator can be written in terms 
of a linear combination of differential operators, and the 
eigenvalues X{p) of adjacent momentum states vary con- 
tinuously, it is intuitive that a distribution function that 
is initially continuous to maintain its continuity as it 
evolves, just like the time evolution governed by a fi- 
nite order differential operator in the usual case. This is 
supported by numerical simulation on the time evolution 
of distribution function. 



IV. CONTINUED FRACTION METHOD 

The tridiagonal nature of the collision operator K 
prompts us to solve the eigenvalue problem by the con- 
tinued fraction method [8, 15]. Hereafter, we omit the 
eigenvector index ji for simplicity of notation. We intro- 
duce the ratio 

F^^iPo)^^^^, i = l,2,3,..., (38) 

which has the advantage of reducing the independent 
c;oniponents in the eigenvalue; equation. In terms of (t>f , 
we need the values of two components of (j)f, to deter- 
mine the values of rest of the other components. But 
in terms of , the value of one component will fix the 
values of the other components. 

We begin by rewriting Eqs. (16a) and (16b) in terms 
of . Since we restrict our consideration to a specific 
Po momentum subspaces, we omit the po dependence on 
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Pf^ for simplicity. In what follows, we will regard as 
a function of A. We get 



1 + n(2 - 2po) 



-A + n(2 - 2po) + n(2 + 2po) 



1 + n(2 + 2po) 

IpoI < 1, 

(39a) 



whereas Eq. (16a') becomes 

^1 (^)--A + n(4)' 

The coefficients are 

l + ?i(2|po;±,| +2) 

n(2|po;±»| - 2) 
1 



Po = ±l. 



x± = 1 



n(2|po;±»|-2) 



„± ^ 1 „ ^± ^ l + ^(2|po:±^|+2) 
' n(2|po;±.|-2) 



(39b) 
(39a') 

(40a) 
(40b) 
(40c) 



When we iterate using Eq. (39b), we get contin- 
ued fraction. In practice, wc truncate the iterations at a 
large enough value of i = 3> 1. This approximation 
reduces the dimensionality of the vector space to 2N + 1 
for pQ ^ 1, and + 1 for = i 1 ■ In this way, the A^-th 
approximant of F^ is 



F±(A) =x±- Ar± + 



X2 - Xr2 



Vn-i 



In the limit of large A'' ^ 1 , we find that 



■''A'' N ^ ^ 



Vn 



_g2;3|po;±jv 



(41) 



(42) 
(43) 



Consequently, Eq. (39b) gives us the boundary value. 



be obtained by iteration, 

1 , 



1 



1 



9 AT 



F, 



N 



Fn 

/ N 



(45) 



This method had been used to solve differential kinetic 
equations [8] and Schrodinger equation [21] too, where 
the continuous nature of the system permits the eigen- 
function expansion in terms of a complete set of orthog- 
onal functions, leading to a tridiagonal structure of the 
equations in terms of the coefficients of expansion. On 
the other hand, due to the discrete nature of the system 
we consider here, we cannot expand the eigenvector in 
terms of continuous orthogonal functions. We instead ap- 
ply the continuous fraction method directly to the com- 
ponents of the eigenvector at discrete set of momentum. 
This method is particularly suitable in solving difference 
kinetic equations with discretely coupled momentum or 
position states. 



V. SPECTRUM 

Due to a reflection symmetry of the collision opera- 
tor under po — Po, we have \^{Pq) = ^)j.{—po), see 
Appendix C. Therefore, we consider only the eigenval- 
ues for po > 0. Fig. 2 illustrates how the eigenvalues are 
determined. The solid curves are F^(X) whereas the dot- 
ted curves are the RHS of Eq. (39a) . The intersections of 
both curves give the eigenvalues of the collision operator. 



A. Accumulation point 




F^CX) 



(l-A)e 



2/3|po;±]v 



(44) 



The continued fraction (41) with the boundary value (44) 
defines F^ as a function of A. By equating the RHS of 
Eqs. (39a) or (39a') to F^ of Eq. (41), we can determine 
the eigenvalues, which are the intersections of these func- 
tions of A. The components of the eigenvector can then 



FIG. 2: For N = 50, T ^ 10 and po = 0. Solid curves 
are F+{X). Dotted curves are the RHS of Eq. (39a). The 
intersections of both curves give the eigenvalues. 

Fig. 3 shows that the spectrum is discrete and bounded 
from below by the steady mode Aq = 0. Furthermore, 
it is bounded from above by a largest value labeled by 
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xlO 



xlOO 



3 - 



1 -■ 



T = 1 



2 r 



1.5 -- 



A 1 



0.5 




T = l 



10 



100 



FIG. 3: Spectrum of the eigenvalues as po varies across < 
po < 1 for three different temperatures, T = 1, 10, 100. The 
top figure is plotted with different vertical scales. The bottom 
figure is a magnified view of the top figure for the interval 
< A < 2. 



A_i below. The most striking feature of the spectrum is 
the existence of a common accumulation point at A = 1, 
which is the special feature in the solution of the differ- 
ence collision operator. 

To see how the accumulation around A = 1 occurs, 
let us consider the matrix K. As mentioned earlier, it 
has a tridiagonal structure according to Eqs. (16a) and 
(16b). For large momentum, or equivalently large i » 1, 
Eq. (16b) approaches 

-Mf- 



i > 1. 



(46) 



The portion of the characteristic equation of K involv- 
ing the large i » 1 elements of the positive branch 0+ 
reduces to (recall the difference of the sign in A from A 
in (17)) 



\K + All = 



A- 1 _ 1 
A- 1 _ 1 
A- 1 



0, (47) 



which gives rise to degenerate value at A=l. The same 
observation occurs for the negative branch <p~ with i <C 
— 1. Therefore, the eigenvalues accumulate around A = 1. 

We recall that we have approximated the continued 
fraction by truncating it at a large momentum po;±Af- As 
we increase TV, new eigenvalues emerge increasingly close 
to the accumulation point, and are distributed evenly on 
both sides of it. On the other hand, the values of the 
smaller and larger values away from A = 1 are not af- 
fected much by the increase of A'', as long as large enough 
iterations N had been carried out. 

In view of the existence of accumulation point and a 
maximum value of A, we label the values in the following 
way. For A < 1, we label them with increasing positive 
integer according to their magnitude, i.e., = Aq < Ai < 
A2 < . . . < 1. For A > 1, we label the values with negative 
integer in decreasing order of magnitude, with the largest 
value of A labeled by A_i, i.e., 1 < . . . < A_2 < A_i < 00. 



B. Band and pairing structure 

Due to the disjointness in the momentum subspaces, 
the eigenvalues in general depend on the representative 
momentum po and vary continuously with it, except the 
steady mode which remains at zero. Consequently, the 
spectrum exhibits a band structure as we vary po, which 
is more prominent for larger values of A. In each momen- 
tum subspace, the spectrum is bounded by the largest 
value A_i (i.e., the smallest value of A), which approaches 
infinity as \po\ — >■ 1- For the boundary values po = ±1, 
we show in Section VIA below that all values with odd 
indices Aj-i, A-1-3, . . ., drop out from the spectrum. A_2 
then becomes the maximum value in this case. The exis- 
tence of a maximum value (i.e., the minimum eigenvalue) 
in this case is distinctively different from the spectrum 
of a differential operator that is usually unbounded from 
above. 

The bottom plot of Fig. 3 also suggests a pairing of the 
values for A < 1. This observation can be explained as 
follows. As a function of A, (41) has infinite numbers 
of poles in the real axis of A. Whenever the variable A 
crosses a pole, the RHS of Eq. (39a) (dotted curves in 
Fig. 2) will intersect (solid curves in Fig. 2) at two 
eigenvalues situated close to the pole on both sides. This 
gives rise to the pairing structure that is more obvious 
for A < 1. 



C. Temperature dependence 

In the low temperature limit T — >■ 0, Eqs. (16a) and 
(16b) go into 



-A0O 



(48a) 
(48b) 
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In this limit, the characteristic equation of K takes the 
form 



\K+~XI\ 



1 A- 1 
1 A _ 1 
A-1 1 



. (49) 



Therefore, there is one eigenvalue at A = 0, whereas the 
eigenvalue corresponding to the accumulation point A = 
1 becomes infinitely degenerate. 

As the temperature increases, the eigenvalues are dis- 
tributed evenly away from the accumulation point. As a 
result, the larger and the smaller eigenvalues are affected 
in an opposite manner by the change in temperature. 
When temperature increases, the magnitude of the larger 
values of A increases whereas the smaller ones decreases. 



VI. EIGENFUNCTIONS 

In the usual eigenvalue problem of a finite order differ- 
ential operator, we consider an equation of the form 



(50) 



where the eigenvalues are independent of p. However, In 
the ID model we consider in the present case, we already 
learned that because of the disjoint nature of the momen- 
tum space, the eigenvalues are generally p-dependent, ex- 
cept for the degenerate steady mode. Consequently, we 
consider an equation of the form 



(51) 



where the eigenvalues now depend on p. As a result, 
the completeness and orthogonality relations only hold 
for components of (j>(p) belonging to the same momen- 
tum subspaces 4>{po;±i)- Hence, in general the idea of an 
eigenfunction with p as a continuous variable does not ex- 
ist in the exact sense of Eq. (50). However, when Xn{po) 
are approximately degenerate, for example for A^ < 1, 
or for A^ > 1 at high temperature, see Fig. 3, one can 
construct approximate eigenfunctions of the collision op- 
erator K. that is at least piecewise continuous. 

In this section, we regard the eigenvalues and eigen- 
vectors as continuous functions of po in the range [—1, 1]. 
In the limit po — )■ ±1, the left side limit (±1~) and the 
right side limit (il"*") of the component of the eigenvec- 
tor at odd value of momentum, podd = ±1,±3,±5, . . ., 
may not be the same. This will be shown to cause the 
disappearance of the eigenvalue at po = ±1 for odd fj,. 
After clarifying this fact, we construct the piecewise con- 
tinuous functions and compare them to the approximate 
eigenfunctions obtained by truncating the collision oper- 
ator to become a second order differential operator at the 
high temperature and large momentum limit derived in 
Ref. [14]. 



A. Disappearance of eigenvalues with odd indices 
when Po = ±1 

An interesting feature of the spectriim occurs in the 
boundary case |po| = 1- In this case, both branches of 
the momentum po;±i coincide, thus reducing the inde- 
pendent components of the eigenvectors by almost half, 
see Fig. l(iii). It is found that the eigenvalues with odd 
indices in the |po| < 1 spectrum, i.e., A±i,A±3,..., dis- 
appear from the spectrum of |po| = 1- As a result, there 
is a discontinuity of (f)^ at podd for these eigenvalues. 

To understand the disappearance, let us consider the 
limit Po — >■ 1 (the other limit po — >■ — 1 can be worked out 
similarly). In Eq. (16a), the factor n(2 — 2po) diverges in 
the limit po — > 1. In order to avoid this divergence, and 
to have Eq. (16a) goes into Eq. (16a') in the same limit, 
we require 



lim 

Po->-l" 



{n(2 - 2po)</>^;o - [1 + n(2 - 2po)](/>+.i} = . 



(52) 



where po — >■ 1 or 1+ means that po approaches 1 
infinitesimally from the left or right side, respectively. 
Since 1 -|- n(2 - 2po) = n(2 - 2po) exp(-|2 - 2po|^), we 
conclude that 



lim 

Po-*-i- 



lim 

Po->-l" 



(53) 



From the notation defined in Eq. (15), Eq. (53) is equiv- 
alent to 



lim 0^(po)= lim (?!)^(-po -h 2) 



= lim J)^{po). 



(54) 



If we consider the limit po — >^ — 1 , we will instead obtain 



_ lim (j)^(po) = _ lim 

po— >■— 1 po— >— 1 



.(-Po - 2) 



lim (pi^ipo) ■ 

po->— 1+ 

(55) 



Through a numerical study, we find that the require- 
ments (54) and (55) are satisfied by solutions that corre- 
spond to eigenvalues with even indices, Aq, A±2, A±4, .... 
These eigenvalues are always finite. As for the solu- 
tions that correspond to eigenvalues with odd indices, 
Ai, A-i-3, . . ., the limit happens to vanish and both re- 
quirements are satisfied. Eq. (16b) then implies that all 
other components of the eigenvectors, (jx^.^, also vanish, 
and this set of eigenvalues drops out from the spectrum 
of Po = ±1. The maximum value is then A_2, which is 
always finite. 

An exceptional case occurs for ji = —1. The limit 
Po ^ 1~ but Po ^ 1 in the components (j)^.fl and (^J.^ 
is non- vanishing, and Eq. (53) is not satisfied. Hence 
the largest value A_i becomes increasing large as po ap- 
proaches 1 from the left side, due to the factor n(2 — 2po). 
An infinite decay rate means that this mode vanishes al- 
most immediately in the time evolution. It is then the 
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second largest value A_2, whieh is always finite, that is 
observed in the time evolution of the state. At exactly 
Po = 1; to avoid the divergence of n(2 — 2pq), the compo- 
nents 0^:0 and have to vanish identically. Therefore, 
A_i drops out from the spectrum at po = ±1- 



Now we show that for fj, 



1,±3,±5, 



has discontinuities in the neighborhood of Podd = 

±1,±3,±5, Using Eqs. (16a) and (16b), we consider 

the sum 

A^(<^m;o + 0^;i) = n{2 + 2po)<^^;0 + n(6 - 2po)</'+i 
- [1 + n(2 + 2po)]</'- 1 - [1 + n(6 - 2po)]<^+2 • (56) 

Since the limit on the both sides of Eq. (53) happens 
to vanish as already shown in the previous paragraph, 
Eq. (56) implies that 

= lim (j)^{p) + _ lim (j)^{p) = Q. 

Po^— 3~ po— > — 3+ 

(57) 

For n = —1^ the limit on both sides of Eq. (53) do not 
vanish, but tlic;y arc; of equal magnitude Init of different 
sign. Hence in the limit po ^ 1~ Eq. (56) again leads to 
Eq. (57). For other poAd, we can then derive inductively 
that 



lim () 



AP) 



lim (t>i^{p) 



(58) 



for all odd /i. In other words, there is a discontimiity 
of (/)p(podd) at Podd for the eigenvalues with odd indices, 
whereas 0^(podd) = 0. This can be seen from the (ii) and 
(iv) plots in Fig. 4. 



B. Piecewise continuous functions 
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FIG. 4: Approximate eigenfunctions of the collision operator, 
only the steady mode (i) is an exact cigcnfunction. Each plot 
is made by a collection of points for momentum subspaces 
connected to po = 0, ±0.1, ±0.2, . . . , ±1, where T = 10. 



constructed in this way. Symmetric (j)^ occurs for 
H = . . . , —3,-1, 0, 2, 4, ... , and anti-symmetric oc- 
curs for /z = . . . , —4, —2, 1,3, 

As discussed in the first part of Section VI, (pfj.{p) be- 
come the eigenfunction of the collision operator when the 
eigenvalues for different momentum subspaces become 
approximately degenerate. As can be seen from Fig. 3, 
the eigenvalues are always close to degenerate for < 1, 

i.e., when /z = 0, 1, 2, As for A^ > 1, the eigenvalues 

are approximately degenerate only in the high tempera- 
ture limit for jj ■ 



-2,-3,-4,.... 



The piecewise continuous functions for each eigenvalue 
indices can be constructed as follows. After normaliz- 
ing 0^;o = (t>^iPa) according to (((/)^(po)|(/)p(jlo)» = 1, 
we are still left with the freedom to choose the sign of 
the representative component 0/t(po)- Once this sign is 
fixed, the signs for the rest of 4>^{po-±i) are automat- 
ically fixed. We use this freedom to set (t>^{po) > for 
Po > 0. Whenpo < O7 the reflection symmetry (Cll) dis- 
cussed in Appendix C gives rise to even or odd functions, 
^^^(^o) = '/>/i(-Po) or -(/)^(-po), respectively, depending 
on the eigenvalue index /x. 

In this way, we obtain even functions (p^, (p) whenever 
H = 0, 2, 4, ... , —5, —3, —1, and odd functions whenever 
H = 1, 3, 5, . . . , — 6, — 4, — 2. It happens that the func- 
tions (f>fi{p) for /i = ±1, ±3, ±5, . . ., can at best be made 
piecewise continuous, due to the discontinuities at Podd, 
which is a consequence of the disappearance of the eigen- 
values from the spectrum for \Pq\ = 1, as discussed in 
Section VIA. In Figs. 4, we give the examples of a few 



C. Comparison with solution of approximate IC 

In this section we will fist quote the results of an ap- 
proximate collision operator obtained in Ref. [14]. The 
approximate collision operator is obtained by truncating 
the series of differential operators of IC at the second or- 
der, supplemented by the conditions T ^ 1 and \p\ ^ 1, 
while maintaining T ^ \p\. Labeling the eigenfunctions 
and eigenvalues as (f> and A^, respectively, the approxi- 
mation results in the following second order differential 
equation 



p dp^ 



dx 

—u 

dp 



2 



(59) 



where the function = Vcq'P^^i in which ipcq = (j>o = 
exp(— p^/2T')/v27rr is the steady state solution. The 
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normalized eigenfunction is 



g-r/4T g-(p-A„T)V4T 



^Jf/2{p/f--X^) 



(60) 



with eigenvalue 



A. = 27^. (61) 

The eigenfunctions satisfy the orthonormahty condition 

/oo 
-00 

A comparison between the spectrum of \^ with the 
solutions of the continuous fraction method in Table 1 
shows that closely approximate the solution obtained 
from the continued fraction method for the representative 
momentum pQ = ±1 for T = 100. The comparison also 
shows that the approximate solution misses the eigenval- 
ues with odd indices in the exact coUision operator. The 
missing modes correspond to piecewise continuous func- 
tions ^odd, which cannot be the solutions of the second 
order differential equation (59). Hence, their omission 
from the spectrum of the approximate continuous colli- 
sion operator is natural. Note that in the approximate 
solutions, the maximum of the non-steady modes of 4>^ 
is shifted from p = 0. Therefore, <f) loses the reflection 

' — M 

symmetry exhibited by along the p = axis. 
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rise to a complete set of biorthogonal basis that belongs 
to the generalized space. 

It is a special feature of the ID nature of the system 
that the momentum and energy conservation permit mo- 
mentum states to be coupled only discretely. This results 
in infinite sets of disjoint momentum subspaces dynami- 
cally independent of one another. The discrete nature of 
the collision operator facilitates the solution of the eigen- 
value problem by the continued fraction method. We 
expect that this method is applicable as well in solving 
other ID models that have discretely coupled space or 
momentum subspaces. 

We find that the spectrum of the collision operator is 
real, discrete and non-negative, and the system mono- 
tonically relaxes towards equilibrium. The most striking 
feature of the spectrum for the difference collision op- 
erator is the existence of an accumulation point. It also 
has an minimimi eigenvalue; in each momentum subspace. 
This is in sharp contrast to the spectrum of the usual 
dissipative differential operators that is unbounded from 
above. 

Except for the steady mode, the eigenvalues of the de- 
caying modes are non-degenerate for different momentum 
subspaces, giving rise to a band structure in the spec- 
trum. This structure originates from the disjoint nature 
of each momentum space. 

In higher dimensions, all the momentum states are con- 
nected and this structure is lost. The general features of 
the results found here show up in other ID models as 
well, such as in the ID quantum Lorentz gas [22]. 

Research is under way to study the time evolution of 
specific distribution functions, with the aim to elucidate 
the consequences of the accumulation point in the spec- 
trum on the time evolution. The accumulation point 
plays a dominant role in the time evolution of the sys- 
tem at low temperature, and it might be related to the 
mechanism of spontaneous emission. 



TABLE I: Comparison of eigenvalues obtained from the con- 
tinuous fraction method \^ for po = ±1 (with missing odd 
indices eigenvalues), with those obtained from the approxi- 
mate collision operator, A^, for T = 100. 



VII. CONCLUSION 

We consider the relaxation process of a quantum 
Brownian oscillator surrounded by a thermal equilibrium 
phonons in a ID chain as a reduced dynamics of the Liou- 
ville operator. The collision operator is a non-hermitian 
operator with dissipative effect that is purely quantum in 
nature. The collision operator for the momentum distri- 
bution function takes the form of a difference operator. 
The eigenvalue problem of the collision operator gives 
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Appendix A: Positivity of eigenvalues 



and po for simplicity) 



Applying the orthogonality relation (30) to Eqs. (16a) 
or (16a'), and (16b), and then using the relations po;±i = 
—po ± 2 and 



1,2,3,..., (Al) 



|PO;±(i+l)l = |PO;±i| +2, 

we can show that the eigenvalues are non-negative. 



-,8Po/2 



n(2-2po)|x+i-X^;o|' 
+ n(2 + 2po)|X;I;i -Xko\ 



+ |2 



+ ^ [e-'^p^./2^2|po;i| + 2)|x+,+i - X+, 
=1 

+ e-^Po;-./2„(2|po._.| + 2)|x-i+i - X-.i 



>0, 



bol < 1 , 



(A2) 



-(1- A)e2^|P"l 



1 ± W 1 



4g-2/3|pjv| 

(1-A)2 



(Bl) 



The vanishing exponential in the square root allows us 
to express the square root as a series of polynomials. As 
a result, the two solutions are 



5^(1) 

■ N 



(1 - A)e 



2mpN\ 



_ _!_ 

^ ~ 1 - A 



(B2) 
(B3) 



F^^ is just the boundary condition for we have made 
use of, see Eq. (44). 

Using Eq. (45), the first solution is 



"at 



:-W\Pn\ 

1-A 



on-1 = ■■■ 



A 



(B4) 



and 



. ^ e-^l>lT^/M2\po■,^^\ + 2)|xJ,+i - Xlf 



>0, 



(A3) 



since the phonon number density n is positive. Therefore, 
the system monotonously evolves towards equilibrium. 
The steady state is a constant column vector, xo\i = 1 
for all i. Hence, it is clear from Eqs. (A2) and (A3) that 
Ao = for the steady mode. 



Appendix B: Solution selected by continued fraction 

In the eigenvalue equation (16a), when any two of (pQ, 
(f)t and (j)^ are chosen independently, the rest of the 
can be determined from the subsequent equations in 
(16b). Therefore, this set of equations should have two 
independent solutions. However, when we write Eq. (16a) 
in terms of the ratio F^ = ^g/'/'f , the continued fraction 
method gives unique F^, and therefore (pf and are 
no longer independent once (po is given. This shows that 
the continued fraction method yields only one of the two 
possible solutions. We will now show that the solution 
that is not produced by the continued fraction method 
can be divergent at large p. Therefore, the method au- 
tomatically discards the unphysical solution. 

For large i = N ^ 1, we can use Eq. (39b) to solve 



for by approximating F^^^ 



F, 



N- 



We then use 



which vanishes rapidly for large N. It is known that 
continued fraction method selects the solution with the 
behavior that as N increases, 4>n either decreases in the 
fastest way or increases in the slowest way [8]. c/)^^^ is 
in fact the solution selected by the continued fraction 
method, and is consistent with the vanishing requirement 
of the eigenfunction at infinity. On the other hand, the 
second solution is 



/,(2) 



(1 - A)0jv- 



For A > 1, (f)j^ can be a large quantity and is not consis- 
tent with the requirement that the eigenfunction should 
vanish at infinity. This solution is not produced by the 
continued fraction method. 

From Eq. (B4), we also learn that if we are looking 
for eigenvalue that is close to A = 1, wc need to iter- 
ate the continued fraction even more times (larger N), 
so that the approximation we used to expand the square 
root in Eq. (Bl) remains valid. On the other hand, one 
may suggest using the ratio F^ = 1/F^ for iteration in 
Eqs. (16a) or (16a'), and (16b). However, it is found that 
solution in terms of F^ produces only the eigenvalues 
with even indices. This is because for the solutions that 
correspond to the eigenvalues with odd indices, <j>o is ei- 
ther exactly zero or close to zero. The ratio F^ = (pf /cpo 
is therefore a divergent or very large quantity. The con- 
tinued fraction method does not generate solutions for 
this situation. 



(1-A) 



(B5) 



Eqs. (42) and (43) to obtain (omitting the ± superscript 



Appendix C: Reflection Symmetry of K 

In this appendix we establish a reflection symmetry 
between the matrices K{Pq) and K{—po), which leads 
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to a relationship between and its mirror reflection 

We start by defining a x -operation tliat acts on a 
square matrix A, and a column vector ip in the following 
ways, 

A^. ^ A_,,_j , ^ V-i • (CI) 

This operation shuffles the components of the matrices, 
but leaves the values of the components intact. 

A closer look at Eq. (16b) shows that it depends only 
on the magnitude of po;i- Furthermore, from Eq. (13), 
we have the relation 

{-po)■,^ = i-in-po - 2i) = -i-lTiPo + 2i) 

= -{PO;-i), (C2) 

which implies 

\{-poh\ = \PO;-i\- (C3) 

Based on these facts, we find that 

(i^x)i,,(-po) = K_i,_j{-Po) = KijiPo) , (C4) 
which is equivalent to 

K''{-p„)^K{po). (C5) 

Now we carry out the x -operation on the eigenvalue 
equation (18) for —pQ. Using Eq. (C5), we obtain 

KiPo) ■ <p; i-Po) = -M-Po)4>; i-Po) . (C6) 



On the other hand, the components of 0^ {—po) is 

?^^;i(-P0) = <t>ij.--i{-po) = <f),i[{-po);-i] = <t>iJ.{-PO;i) , 

(C7) 

where we have used Eq. (C2) in the last equality. In 
component's form, Eq. (C6) becomes 

Kij{po)(t)^,{-po■j) = ->^ij.{-po)(t>ix{-Po;i) ■ (C8) 

Therefore, (f>f^{—pn-i) is the eigenvector of K with eigen- 
value Ap(— Po)- Comparing this with the component's 
form of the right eigenvalue problem (18), 

Kij{po)(t)n{PO;j) = -V(Po)0/:i(PO;i) , (C9) 

and using the fact that the eigenvalues are real, non- 
negative and non-degenerate, we can match the eigenval- 
ues according to 

V(-Po) = V(Po) . (CIO) 

Furthermore, 4>fj,{—po:i) is related to 4>fj,{pa-i) up to a 
phase. Since K and the eigenvectors are real, the phase 
can only take the real value ±1. We then conclude that 

<Pf^{-pO;i) = ±0/x(PO;i) • (Cll) 

This result enables us to construct a function 4>^{p) in 
Section VI B that is piecewise continuous in p, and is ei- 
ther symmetric or anti-symmetric with respect to mirror 
reflection in p. 
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